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SUBROUTINE  IDENTIFICATION 

A.  Title:  Compute  Satellite  Position  in  Cartesian  Coordinates  (CARTES) 
Ident.  084,  Mod.  EA. 

B.  Programmed:  D.  Leder,  EI&P,  MSVD  General  Electric  Company. 

C.  Documented:  15  October  1962,  J.  LaVine,  System  Development 
Corporation. 

PURPOSE 

CARTES  coniputcs  the  position  of  a  satellite  in  Cartesian  Coordinates 
using  given  orbital  position  parameters, 

USAGE 


Calling  Sequence. 
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1*12 

NOP 

ZRO 

B 

L(u>  ) 

P 

L+13 

NOP 

ZRO 

B 

I*(r) 

l*lk 

NOP 

ZRO 

B 

L(x) 

L+15 

NOP 

ZRO 

B 

L(y) 

L+l6 

NOP 

ZRO 

B 

L(z) 

L+17 

NOP 

ZRO 

B 

L(5c) 

L+18 

NOP 

ZRO 

£ 

L(y) 

L+19 

NOP 

ZRO 

B 

L(i) 

L+20 

NORMAL  RETURN 

The  address  protion  of  each  lower  instruction  refers  to  the 
location  in  core  of  the  particular  parameter.  The  B  term  may 
he  zero,  or  any  legal  index  register  designation. 

B.  Input  Parameters 

The  following  parameters  must  he  in  the  locations  specified  by  the 
calling  sequence.  All  values  are  to  be  in  floating  point  format. 


PARAMETER 

UNITS 

DESCRIPTION 

J 

N.D. 

Coefficient  of  second  spherical 
harmonic  in  Earth's  gravitation  field, 

R€ 

feet 

Earth's  equatorial  radius. 

g 

feet/sec2 

Acceleration  due  to  gravity. 

% 

radians/sec. 

Earth’s  rotation  about  polar  axis. 

i 

radians 

Inclination  of  orbital  plane  to 
equatorial  plane. 

l 

0 

feet 

Semi-latus  rectum  of  mean  orbit  at 
time  of  first  ascending  node. 
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PARAMETER 

4 

*  03 


*  t 

1 

*  8. 

*  e 

*  03 

] 

*  r 


UNITS  DESCRIPTION 

radians  Longitude  of  first  ascending  node. 

radians  Total  angular  travel  from  first  ascending 
node. 

seconds  Time  measured  from  nominal  time  of  launch. 

feet  Semi -major  axis  of  mean  orbit. 

N.D.  Eccentricity  of  mean  orbit, 

radians  Argument  of  perigee  of  mean  orbit 

feet  Geocentric  distance 


*  These  position  parameters  are  obtained  from  the  output 
of  either  the  KH1  or  KH2  subroutine. 

C.  Results 

The  outputs,  which  consist  of  the  components  of  the  radius  and 
velocity  vectors,  are  stored  as  specified  by  the  last  six  address 
in  the  calling  sequence.  These  values  are  in  tee  form  of  Cartesian 
coordinates  which  represent  a  geocentric  inertial  system  with  the 
positive  Z  axis  pointing  North  along  the  polar  axis,  the  positive  X 
axis  along  the  right  ascension  of  the  Greenwich  meridian  at  nominal 
time  of  launch,  and  the  Y  axis  90°  East  of  the  X  axis.  They  are 
defined  as  follows: 


PARAMETER 
x,  y,  2. 

y,  '> 


UNITS  DESCRIPTION 

feet  Coordinates  of  satellite  position  (floating 

point) 

feet/sec.  Velocity  components  in  the  positive  ,  y,  z 
direction  (floating  point) 


RESTRICTIONS 

A.  Accuracy: 

The  accuracy  of  this  program  is  affected  only  by  the  limitations  of 
floating  point  arithmetic  and  the  results  obtained  from  the 
referenced  subroutines.  This  assumes  the  input  orbital  elements 
are  within  a  valid  range, 
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B.  Environment: 

CARTES  used  the  subroutines  SIN,  COS,  SQRT,  and  SUBERR.  The  program 
uses  and  restores  index  registers  1  and  2.  No  Reference  Fool  items 
are  used. 

C .  Error  Checks : 

CARTES  makes  no  validity  checks  on  inputs  from  the  user  program  or 
values  obtained  from  other  subroutines.  The  only  type  of  error 
exit  provided  is  to  the  SUBERR  subroutine  upon  return  from  the  SIN, 
COS,  or  SQRT  subroutines. 

TIMING 


The  time  required  for  CARTES  to  cycle  is  3972.2  microseconds  plus  the 
time  required  for  the  associated  subroutines.  This  figure  disregards 
the  possibility  of  a  premature  exit  caused  by  an  error. 

STROAGE  REQUIREMENTS 


A  Space  Allocation: 
Program 
Constants 
Temporaty  Storage 
TOTAL 


121  cells 
6  cells 
4l  cells 
'  168  cells 
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